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SUMMARY 

A review  is  given  of  analytical  methods  for  predicting  the 
muzzle  velocity  of  compressed  gas  guns.  Steady  adiabatic 
expansion  theory,  quasi-steady  theory  and  unsteady  theory  are 
all  examined,  with  emphasis  placed  on  simplified  analytical 
methods  applicable  to  subsonic  and  transonic  muzzle  velocities. 
New  formulae  are  obtained  and  compared  with  numerical  results 
from  a full  theoretical  treatment  of  the  gas  flow.  The 
comparisons  show  that  the  new  analytic  results  are  remarkably 
effective  and  can  therefore  be  used  for  gun  design.  The 
effects  of  atmospheric  counter  pressure  are  treated  separately 
and  a new  method  of  correcting  the  predicted  velocity  is  given. 
An  analysis  of  data  from  actual  subsonic  guns  highlights  the 
importance  of  muzzle  pressure  reflections  in  reducing  the 
retarding  pressure  at  low  subsonic  speeds.  The  observed 
muzzle  velocities  are  about  10%  below  the  muzzle  velocities 
predicted  by  theory.  Finally,  nitrogen,  helium  and  hydrogen 
are  compared  as  potential  driver  gases  for  transonic  guns. 
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1.  INTRODUCTION 

This  report  is  concerned  primarily  with  the  gas  dvnainic  theory  underlying  gas 
guns  operating  with  muzzle  velocities  in  the  transonic  range.  A good  reference 
to  the  fundamental  theory  of  unsteady  gas  motions  is  the  review  of  Seigel (ref . 1) 
who  has  applied  the  theory  to  two-stage  h>qpervelocity  launchers  in  which  a moving 
piston  is  used  as  a gas  compressor.  More  than  this,  Seigel  has  given  a compre- 
hensive coverage  of  the  theory  of  single-stage  gas  guns  such  as  are  considered  in 
this  report.  The  hulk  of  tlie  design  cliarts  in  reference  1 apply  to  single-stage 
guns  and  are  therefore  relevant  to  the  present  investigation.  A few  of  the  many 
papers  concerned  with  gas  dynamics,  hypervelocity  launchers  and  the  use  of  moving 
pistons  as  gas  compressors  are  given  in  references  (2, 3, 4, 5, 6, 7') . Reference  2 
is  included  as  an  example  of  a book  providing  an  introduction  to  unsteady  gas 
dynamics.  References  3,4  and  S contain  additional  bibliographies.  Reference  6, 
on  the  other  hand,  is  concemed  solely  with  the  use  of  a moving  piston  to  compress 
gas  to  a very  higli  pressure.  Reference  7 describes  the  application  of  a piston 
comjiressor  to  drive  a hN'personic  wind  tunnel. 

Interest  in  gas  guns  began  at  Weapons  Research  Establishment  ten  years  ago. 
During  an  investigation  into  tlie  flight  of  lifting  projectile  shapes,  it  was 
suggested  that  a gas  gun  could  fire  single  projectiles,  mounted  in  sabots,  and 
that  measurements  could  be  made  of  their  subsequent  flight.  A gun  of  127  mm 
bore  and  capable  of  operating  at  pressures  up  to  3000  kPa  was  constructed  and 
has  been  in  almost  continuous  use  since  its  first  firing  late  in  1967(rcf.8). 

Sabot  speed  is  measured  at  the  muzzle  of  the  gun.  The  gun  has  supported  tasks 
ranging  from  the  statistical  study  of  projectile  impact  patterns  to  assisting 
in  the  development  of  a W.R.E.  image  svneh roni zat i on  camera.  Early  in  1973, 
a larger  gas  gun  of  384  mm  bore  was  commissioned  to  study  the  flight  of  larger 
projectiles  and  clusters  of  small  project! les (ref . 9) . Tin  s gun  operates  at 
pressures  up  to  1000  kPa.  The  two  gas  guns,  together  with  associated  instru- 
mentation, provide  a facility  for  aerod\aiamic  research(ref.  10,11)  £md  the  develop- 
ment of  small  munitions  in  a c lose-to-the- laboratory  environment  where  remote 
and  extensive  equipments  such  as  those  at  Woomcra  in  South  Australia  are  not 
necessary. 

The  gas  guns,  a control  and  recording  building  and  pressurizing  equipment  are 
located  at  one  end  of  an  area  approximately  1100  m by  600  m in  size.  Tlie  ground 
instrumentation  includes  some  fixeil  wide-angle  ballistic  cameras  which  are  suit- 
able for  trials  at  night  to  determine  projectile  position,  attitude  and  velocity 
throughout  flight.  Tivo  lights  flasliing  at  a frequency  of  about  50  Hz  are 
carried  in  the  projectile  and  the  light  flashes  recorded  on  the  ballistic  camera 
plates.  During  day  firings  performance  of  the  sabot  and  subsequent  flight  of 
the  projectile  may  be  recorded  by  high  speed  cameras. 

The  W.R.E.  gas  guns  continue  to  bo  used  extensively.  For  this  reason,  pre- 
liminary consideration  has  been  given  to  the  design  of  a gun  capable  of  much 
higher  speeds.  Since  acrodviiamic  I'roblems  of  a varied  nature  frequently  occur 
at  transonic  speeds,  ;iny  such  gun  should  be  capable  of  testing  projectiles  at 
speeds  up  to  at  least  a Mach  inunber  of  1.2,  and  preferably  up  to  a Mach  number 
of  1.5  or  so.  Tlie  purpose  of  tlie  present  report  is  to  provide  the  necessary 
theoretical  design  b;isis.  As  will  be  seen,  the  report  relics  heavily  on  work 
carried  out  by  other  investigators. 

We  begin  with  the  simple  flow  model  described  by  Perfect (ref . 12) . Tliis  model 
is  extended  to  higlier  subsonic  speeds  in  Section  2.  For  supersonic  gas  speeds, 
a different  model  is  necessary  and  one  based  on  unsteady  expansion  waves  is 
described  in  Section  3.  Throughout,  the  emphasis  in  on  simplified  analytical 
methods  which  may  be  applied  readily  by  gun  designers.  Section  4 introduces 
the  retarding  effect  of  atmospheric  counter  pressure  acting  on  the  front  face  of 
a projectile.  Tlie  perform;uice  of  existing  subsonic  gas  guns  is  examined  in 
detail  in  Section  5 in  order  to  assess  the  precision  with  which  muzzle  velocities 
can  be  predicted.  A summary  of  the  recommended  design  procedures  is  given  in 
Section  6.  An  example  is  presented  in  Section  7,  to  compare  hydrogen,  helium  and 
nitrogen  as  driver  gases  for  a transonic  gas  gun.  Finally,  conclusions  are  given 
in  Section  8. 
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The  reader  interested  solely  in  design  application  should  begin  by  reading 
Section  1.1,  followed  by  Sections  6 cuid  7. 

1. 1 Steady  internal  flow  model 

Notation  for  the  analysis  of  gas  gun  performance  is  illustrated  in  figure  1. 
A projectile  of  mass  m is  accelerated  by  compressed  gas  along  a barrel  of 
length  Li . Initially,  the  gas  pressure  is  applied  to  the  projectile  by 
opening  a quick-acting  valve  or  bursting  a diaphragm.  The  equation  of 
motion  of  the  projectile  is 


m 


d^x 

dt^ 


(Pj  - Pj.)-  Ai 


(1) 


wliere  p^  is  the  accelerating  or  driving  pressure  and  p^  is  the  retarding 
pressure. 


Gas  conditions  before  firing 

Reservoir  Barrel 


Pressure 

Po 

Pi 

Speed  of  sound 

a I 

Ratio  of  specific  heats 

To 

7i 

L— ... I..  

RESERVOIR 

Ho 

L,  .-| 

Volume  Vo 

Gas  mass  G 

BARREL 



PROJECTILE 

Mass  m 

Velocity  u 


1 


D, 


Reservoir  area  = (7r/4jDo 
Reservoir  volume  Vq  = AqI^ 
Barrel  area  A,  = (7r/4)D|^ 

Volume  ratio  K = Vq/A,  L| 
Muzzle  velocity  U = (u) 


Mass  ratio  G/m 


Po^o  _ 
m 


7o.k. 

mao 


Figure  1 . 


Notation 
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Perfect (ref. IJ)  presents  a simple  solution  to  equation  (1)  by  using  a 
steady  flow  model.  The  driving  pressure  p^  is  calculated  from  an  adiabatic 

expansion  of  the  reservoir  gas  from  its  initial  volume  Vq  to  the  volume 
(Vq  + Ai  xl . At  tlie  same  time,  the  retarding  jjressure  is  taken  to  be  the 
same  as  the  initial  pressure  in  the  barrel.  Therefore 


Pd  = Po 


Vo /(Vo 


+ Aj  x) 


7o 


(2a) 


and 


]\  = Pi. 

Hie  solution  to  equation  (1)  is  now  found  to  be 


(2b) 


2poVo 
(7o -l)m 


1 + -21 
Vo 


A,  L,V 


(7o-l) 


2pi  A,  L, 


m 


or 


2]>oAi 

K 

1 

1 

c 

1 

+ 

f 

m 1 

1 

_ C.7-T, 

1 

^ ' 

7^ 

2pi  Ai  U 

m 


(3) 


where  k = (Vo/AiLi)»  the  ratio  of  the  reservoir  volume  to  ttie  barrel  volume. 

Tlie  analysis  to  follow  will  begin  by  ignoring  the  effect  of  the  retarding 
pressure  since,  tor  transonic  guns,  the  initial  reservoir  pressure  is 
expected  to  be  very  much  greater  than  the  initial  retarding  pressure.  With 
this  in  mind,  it  is  convenient  to  write  equation  (S)  in  the  alternative  form 


2 Pi  A|  L| 

0 ■ 

III 


(V) 


Equation  (4)  will  be  used  in  subsequent  analysis  of  atmospheric  counter 
pressure  effects. 


2.  QUA.SI-.STliADY  I.VTERNAL  l-l-OW  MODEL 

llie  steady  flow  model  of  Section  1.1  above  relics  on  an  adiabatic  expansion 
of  the  reservoir  gas.  lliis  means  that  tlie  Mach  number  of  the  flowing  gas  is 
always  assumed  to  be  small.  It  is,  howev'er,  stra i ght fonvard  to  derive  a 
correction  to  eejuation  (2a)  and  so  develop  a quasi-steady  model.  Such  a model 
is  discussed  in  reference  9. 
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2.1  General  tlieory 

Following  reference  9,  the  basic  idea  is  to  regard  equation  (2a)  as  giving 
the  stagnation  pressure.  A new  ex]iression  for  the  driving  pressure  p^  is  then 

found  using  st;indard  compressible  flow  formulae  such  as  are  given  in  reference 
2 for  example.  We  find  that 


I’d  = Ps 


70/(70-1) 


(5a) 


whe  re 


'■s 


Po 


I'o  / (I'o  Ai  xj 


7o 


(5b) 


;md  a^,  the  stagnation  speed  of  soujul,  is  given  by  the  square  root  of  the 
temperature  ratio  found  from  the  adiabatic  expansion.  i.e.  (a^/ao )= (T^/T® ) * 


ai'id  .;o 


ao 


i 0 / (Vo  + A|  X I 


(7o-l) 


(Sc) 


Eqracion  il),  with  the  retarding  iircssure  p omitted,  now  leads  to  the 
eejuation  of  motion  ^ 


mu 


du 

dx 


A,Pd  . 


(8) 


rh  1 ■;  equation  c.'uinot  bo  solved  anal  yt  i ca  1 1)'  as  the  variables  u and  x arc  not 
separated,  because  of  tlie  fcnait  of  ecpiation  (Sa) . Numerical  solution  on  a 

computer  is  sim])le  ;md  straightforward.  .-in  malytic  solution  can,  iiowever, 
be  foimd  by  seeking  a low  Mach  niunber  solution.  llten , from  cipiation  (5a) 
we  wri'e 


I’d  = Ps 


7o 

1 


Pj  = po  I !()  / (Vo  + A,  x) 


7o 


7oPoU  j y , (y 

■ '0  / 1^0 

2;i^  I 


+ A,  x) 


using  equations  (5b)  and  (5c). 
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in  agreement  with  equation  (oi  wlien  pi  •■=  0 and  fpo Ai  Li /ma5 )'  is  small. 
Denoting  the  muzzle  velocity  from  equation  (31  by  Uy , we  have 


(Uo/aoV  = 2 (po  Ai  Li /ma5  ) . 


Hence,  when  V'o  is  very  large,  the  result  can  be  written  as 

= y'7o)  ■ . ' - ^ (Do/;io  1; 

or,  since  wo  have  assumed  tiiat  the  Macli  number  is  small, 

u'/a,^  = (U'/.a;,l  . 

Hence  we  obtain 

U/ao  = (Uo  / ao  1 .i 

I 


1 - ^ . 1115/^1^1 


1 - . (Ho/a^l  I . (81 


This  result  shows  in  a very  simple  way  the  magnitude  of  the  quasi-steady 
correction  in  terms  of  the  muzzle  velocity  Ho  calculated  from  the  steady 
flow  model.  In  the  general  case,  equation  (31  shows  that  this  is 


1 

2 poAi  Li 

1 

2 

K - 

1 

z 

Uo  / = j 

_ m.io  _ 

..  To-l_ 

L J 

(91 


WRE-TR-1877(W) 


! 


6 - 


Equation  (8)  is  compared  with  the  exact  result  from  equation  (7)  in  table  1. 

The  performance  of  ecfuation  (8)  is  surprisingly  good  even  at  high  subsonic 
speeds.  The  effect  of  the  parameter  K is  seen  to  be  of  secondary  importiincc. 
Consequently,  equation  (8)  is  taken  to  lie  a satisfactory  first  approximation 
for  all  values  of  K and  y„ . Equation  18)  is  recommended  for  use  in  design 
calculations. 

/\n  alternative  method  of  allowing  for  quasi-steady  effects  is  to  use  a 
simplitied  fonii  of  the  "Pidduck-Kent  .S]iecial  .Solution"  discussed  in  Appendix  1. 


FABTL  1.  QUASI -.STEADY  MAtTl  NUMBER  CORRECTION 


Uncorrected 
mach  number 
(Uo/aoi,  from 
equation  (9). 

Ratio  of 
specific 
heats,  7o . 

Corrected  macli 

number  (U/ao) 

Prom  equation  (7) 

From 

K=1 

K=2 

K=oo 

Mean 

equation  (8) 

1.1 

0.483 

0.482 

0.483 

0.483 

0.483 

0.5 

1.4 

0.475 

0.472 

0.479 

0.475 

0.478 

5/3 

0.467 

0.461 

0.475 

0.468 

0.474 

1.1 

0.693 

0.692 

0.696 

0.694 

0.692 

0.75 

1.4 

0 . 669 

0.662 

0.682 

0.671 

0.676 

5/3 

0.645 

0.6  29 

0.670 

0.648 

0.662 

1.  1 

0.872 

0.869 

0.877 

0.873 

0.  863 

1.0 

1.4 

0.822 

0.805 

0.848 

0.  825 

0.825 

.5/3 

0.  773 

0. 740 

0.824 

0.779 

0.792 

2.2  Barrel  density  effect 

An  additional  effect  may  be  incoriio rated  in  computer  versions  of  the  quasi- 
steady model.  By  allowing  for  the  Jo<:i'ease  in  gas  density  as  Mach  number 
increases  a cliange  in  the  mass  of  gas,  ;ind  therefore  the  pressure,  in  the 
reservoir  can  be  calculated.  Tlie  starting  point  for  our  analysis  is  the 
compressible  flow  density  re  lat  i onsliip 


u' 

s 


l/(7o-l) 


(10a) 


We  now  assume  that  the  gas  in  the  reservoir  is  at  rest  since,  for  gas  guns, 
the  reservoir  arc<i  Ad  is  considerably  greater  tluui  the  barrel  area  Ai  . Ttie 
gas  in  the  barrel  is  taken  to  be  moving  with  the  projectile  velocity  u. 

I'he  mass  of  gas  in  the  barrel,  when  tlic  projectile  has  travelled  a distance 
X,  is  pA|  X which  may  be  written  as  p^A| [ x (p/p^) I . Thus,  in  calculating  the 

steady-state  adiabatic  expansion,  tlie  reservoir  liehaves  as  though  the  pro- 
jectile had  only  moved  a distance  px/p^.  i.e.  the  stagnation  pressure  in 

the  reservoir  is  given  by 


" Po 


Vo/j  V„ 


4- 


Ai  x(p/p^)j 


7o 


(10b) 
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;md  the  speed  of  sound  is  j;u'eii  by 


iio  . 


V() 4 V'o  + A|X(p/p^)  j 


(7(1  - 1 ) 


(10c) 


liquations  (10b)  imd  (10c)  replace  equations  (sli)  ;uid  (Sc).  l;<)uation  (5a) 
is  no\'i  combined  witli  equations  (10)  to  j;ivt'  the  driving  pressure  required 

in  the  eiiuatii'n  of  motion  ((>). 

Some  computer  calculations  have  been  made  (or  actual  yas  K'U''’-  ^I'o 
results  show  tluit  the  density  efloct  gives  rise  to  relativelv  small  increase’s 
in  muzzle  velocity.  'I'ltc  effect  is,  of  course,  more  pronounced  as  the  gas 
Mach  number  increases.  licpiation  (8)  tliereforc  continues  to  provide  a satis- 
factor_\-  analytic  approximation  for  (iredictions  of  muzzle  velocity. 


3.  IJ.MSTliADV  nXPA.NSION  WAVli  MODE!. 


As  the  gas  Macli  number  increases  to  sujiersonic  speeds,  unstead)’  flow  effects 
become  more  and  more  important.  Put  another  wa>',  tlio  time  taken  for  pressure 
waves  to  travel  from  the  barrel  to  the  end  of  tl\e  reservoir,  and  back  again, 
is  no  longer  small  compared  with  the  time  taken  for  the  projectile  to  move  down 
;md  out  of  the  barrel.  The  quasi-steady  model  of  section  2 impilicitly  assumes 
that  a lot  of  wave  interaction  takes  place.  By  this  means  a nearly  steady 
flow  j)attcrn  is  set  up.  Iv'e  now  turn  to  a situation  where  wave  interactions 
arc  not  import;u\t  and  the  flow  field  is  dominated  by  the  simple  unsteady 
expansion  wave  set  up  by  the  moving  projectile. 

3.  1 fonstiuit  area  c;ise 


The  constant  area  case,  Ao  - Ai  , is  considered  first.  lltcrefore,  the 
standard  e()uations  of  one-dimensional  gas  dvaiamics  ’apply.  In  particular, 
the  pressure  behaviour  in  a centred  cxpiansion  wave,  often  referred  to  as  a 
simjile  wave,  is  of  direct  use.  Keferences  1 ;uid  2,  for  example,  provide 
further  information.  The  formula  relating  pressure  and  velocity  in  an 
un s t e a d y e xp ;ui s i on  is 


p/po 


27o/(7o-1) 


(11) 


I'his  formula  is  exact  up  to  the  time  of  arrival  of  the  first  pressure 
reflections  from  the  back  of  the  reservoir.  liquation  (11)  shows  that  the 
maximum  possible  speed,  or  tite  escape  speed,  is  given  by  1 - ((7o -l)u/2ao ) = 0, 
i.c.  the  maximum  possible  speed  is  2ao/(7o-l).  Clearly,  gas  speeds  arc  not 
restricted  to  a maximum  of  Mac)i  1.  Tor  7o=  1.4  for  example,  tlie  maximum 
possible  speed  based  on  the  speed  of  sound  ao  is  Mach  5. 

Equation  (11)  is  well  known  in  analyses  of  shock  tube  performance.  It 
describes  the  initial  reduction  in  pressure  in  the  high  pressure  section. 

Wl\en  coupled  with  a shock  wave  in  the  low  pressure  section  a complete 
description  of  the  initial  flow  field  is  available.  The  motion  of  the 
contact  surface  separating  the  gas  in  tl^e  high  pressure  section  from  the 
gas  in  the  low  pressure  section  is,  in  some  ways,  analagous  to  the  motion 
of  a projectile  in  a gas  gun.  E(iuation  (II)  describes  the  pressure 
accelerating  tl>c  contact  surface  while  the  corresponding  shock  wave  formula 
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Oi'scr tlifs  the  retarding  pressure.  i'lie  major  difference  is  sim|)ly  that  the 
contact  surface  has  zero  mass  (m=0)  and  is  tlierefore  accelerated  u|)  to  maximum 
speed  i nsttuitaneous ly.  The  present  analysis,  in  contrast,  is  concerned  with 
predicting  the  velocity/distance  history  of  projectiles  with  non-zero  mass. 

Returning  now  to  equation  CD),  this  can  be  coupled  with  the  equation  of 
motion  (6)  to  give 


mu 


du 

dx 


= Ai  j)o 


u_ 

ao 


27o/(7o-l) 


whence  the  muzzle  velocity  U is  found  from 


poAi  I..1 

mao 


2 

J 

^ (IJ/ao)  - (2/(7o  + 1))  I 

7o+l 

1 

111  - ((7„-l)U/2ao)] 

(12) 


llie  formula  suffers  from  the  disadvantage  that  U is  given  implicitly.  for 
design  purposes,  tut  explicit  formula  showing  the  variation  of  U with  the 
driving  partuneter  poAiL|/mao  is  more  convenient.  Such  a formula  is  investi- 
gated in  the  next  section. 

3.2  Logarithmic  approximation  for  muzzle  velocity 

An  empirical  fonnula  giving  IJ  explicitly  at  transonic  and  supersonic  speeds 
is  now  sought.  Tlie  beliaviour  of  equation  (12)  over  the  range  U/ao  = 1 to 
U/ao  = 2 immediately  suggests  a logarithmic  or  nearly  logarithmic  dependence 
of  1)  on  po A|  1.^ /mao  . A number  of  logarithmic  formulae  were  investigated. 

The  simplest  form  was  found  to  be 

U/ao  = In  I C . poA|Li/m;i^| 

-70 

witli  C pro^iort  ional  to  7o  • Die  performance  of  this  formula  for  ('.  = S7o 
and  t;  = 67o  is  shown  in  table  2.  The  general  performance  is  surprisingly 
good  both  at  supersonic  speeds  and  over  a very  wide  range  of  values  of  7o • 
Overall,  C = 67()  seems  to  be  the  best  choice.  The  recommended  design 
fonnula  is  therefore 

U/ao  = In  I 67oroAi  Li /mao  I • (13) 

27o 


L 
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TABU;  2.  UNSTEADY  EXPANSION  FORMULAE 


Ratio  of 
specific 
heats 

7o 

r ■■ 

Mach 
number 
U/ ao 

PoA|  L| 

mao 

(from  eq. (12) ) 

Simplified  formula  for  mach  number 

( 2^)  I’o^' ^ 

Mach 

number 

U/ao 

. . 

c = syl 

c = eyh 

0.5 

o 

00 

0.05 

0.13 

0.5 

0.75 

0.501 

0.50 

0.59 

0.75 

1.1 

1.0 

1.09 

0.86 

0.94 

1.0 

1.5 

3.  77 

1.42 

1.50 

1.5 

2.0 

o 

07 

04 

1.89 

1.97 

2.0  1 

0.  5 

0.206 

0.25 

0. 32 

0.5 

0.  75 

0.612 

0 . 64 

0.71 

0.75 

1.4 

1.0 

1.47 

0.95 

1.02 

1.0 

1 . 5 

6.50 

1.48 

1 . 55 

1.5 

J.-'V, 

25.84 

1.9S 

2.04 

2.0 

0.5 

0..35 

0.41 

0.5 

0.  75 

0.  750 

0.70 

0.76 

0.75 

5/3 

1 0 

2.02 

1.00 

1.05 

1.0 

1 5 

12.75 

1.55 

1.61 

1.5 

20  j 102.0 

2.  18 

1 

1 

1 

! 

I 

1 

1 

1 

2.0 

3.3  Area  cEaiig  ; at  barrel /rose rvo i r Junction 

The  queition  now  is  wheliier  equation  (ISJ  can  be  generalized  to  the  usual 
gas  gun  siaiation  where  Ao  > Ai  ;ind  gas  velocities  in  the  reservoir  are  much 
smaller  thin  in  the  barrel.  Fortunately,  the  barrel  entry  sonic  approxi- 
mation, as  described  in  rcFeiencc  1 for  example,  provides  a basis  for  an 
analytic  approximation.  Because  the  local  gas  velocity  at  the  beginning  of 
the  barrel  afiproaches  sonic  speed  witli  increasing  time,  the  fundamental 
assumption  made  in  calculating  the  muzzle  velocity  is  that  the  flow  is  always 
sonic  at  the  barrel  entry.  In  order  to  show  how  the  formula  has  been 
derived  we  begin  by  examining  the  infinite  area  reservoir  in  which  the  gas 
is  at  rest. 

Since  tlic  unsteady  expansion  flow  we  are  considering  is  a simple  wave, 

the  Ricmann  variable  a + ('yo-l)u/2  is  const;int  throughout  the  wave.  At  the 

barrel  entry,  we  take  u = a = (a)  . ;md  the  Riemann  variable  has  the  value 

' ' ^sonic  

((7o‘^l)/2)  (uj  . which  can  he  expressed  as  ao  \/fT'o  + lJ/2  since  (a)  . = 

_ _ sonic  ' V V u / V 'sonic 

;io  \^'(To  > 1 1 is  an  isentropic  flow.  Hence 


a + (l/2)(7o-nu 


ao\/  (7o  + l)/2 
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Rut,  in  isontropic  flow,  p/po  = (ii/.u, ) ^ ^ ^ so  that 


p/po 


'/(7o>T)/2  - C(7o-l)ii/2£k,) 


27o/(7o-1) 


Writing  tlte  tenn  \'^(7o+TV2  as  1 + ( v^C7()+lV2  - 1)  shows  that  the  term 

(.v^(7o  + l)/2  - n is  tlic  sole  addition  to  equation  (11)  resulting  from  the  use 

ot  tile  sonic  entry  approximation.  Following  Sc  i ge  1 ( ref . 1)  an  inverse  area 
ratio  interpolation  is  introduced  in  tliis  additional  term  to  account  for 
arbitrary  values  of  A, /A,,  . The  pressure  result  is 


p/po 


1 . z - ) - 

A.  / c 


27o/(7o-1) 


where 


)1  - (A,/Ao  )i  . i v/{7o^-i)A2  - Ij  . 


Comparison  with  ecjuation  (11)  shows  that  the  counterpart  of  equation  (13)  can 
now  be  found  by  replacing 


Po  by  p„  I 1 + z| 


27o/(7o-l) 


and  ao  by  ao  [ 1 + Z] . Hence  the  generalized  fonn  of  equation  (13)  becomes 


In  [ 1+Z  1 . 


This  result  can  be  simplified  by  t ak i ng  \/(7o  + 1 ) /2  to  be  equal  to  1+ ((7o - 1) /4) 

wliich  strictly  apjilies  when  (7o  ~ 1)  is  "small",  and  expanding  the  logarithm 
in  a similar  way.  i.e.  we  take 

In  [ 1+Z  ] = J (7o-l)/4i  . ) 1 - (A,/Ao)J  . 


Omitting  terms  multiplied  by  (7o-l^,  the  final  result  is  therefore 


11/ ao 


J_ 

47o 


1 - 


2^ 

Ao 


27o 


In 


b7oPoA|  Li  /m;io 


(14) 


ixpiation  (14)  is  the  basic  design  formula  that  was  being  souglit.  In  the  next 
section,  it  will  be  showm  that  tliis  formula  works  surprisingly  well. 

lor  future  reference,  we  note  that  the  driving  pressure  counterpart  to 
equation  (14)  is 


L _ 


11 
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P/Po 


I 


>0-1 

4 ! 


1 


' lajill  iL 

2 / ao 


27o/(7o-1) 


(151 


An  alternative  pressure  t'onnula  clue  to  Seigel  (ref . 11  is  given  in  Appendix  II. 

3.4  Comparison  with  exact  numerical  results 

Seigel(ref . 1)  has  given  a large  number  of  design  charts  based  on  numerical 
integration  of  the  gas  dviKimic  differential  equations  that  describe  the 
unsteady  gas  Flow.  These  charts  provide  exact  results  which  can  be  used 
to  assess  the  usefulness  of  ecpiation  (141.  Figures  2,  3,  4 present  results 
for  7u  = 1.1,  1.4,  5/.3.  llie  broken  straight  lines  are  the  results  given 

by  equation  (141.  Fhe  results  for  Ao/A|  = °°  can  be  compared  directly 
with  Seigel's  results  for  Ao/A|  - 2.5  since,  from  equation  (141,  the  difference 
between  the  Ao/A|  = °°  case  and  the  Ao/A|  ■=  25  case  must  be  extremely  small. 


n 

■To 


!.5 


2.0L 


/V)/A, 

0 25.0 

X 4.0 

o 1.0 

Data  from  Seigel , 
figures  20, Ref. 1 


Theory  A<, /A,  = °° 


I .5 


0/  / 


1.0 


0.5 


ay  ^ 

/// 


Theory  A^/A,  = 4.0 

Theory  Ao/A|=  1.0 
(equation  (1411 


(TTiT 


1.0  2.0 


5.0 


10.0 


/'  Po  A|  L| 

\ mafi 

Figure  2.  Theory  for  infinite  reservoir  gun,  7o  = 1.1 

For  7o  =1.1,  figure  2 compares  equation  (141  with  Seigel's  exact  results. 
The  exact  results  apply  to  the  case  where  the  reservoir  is  effectively 
infinite  in  length  l.e.  it  is  long  enough  for  wave  reflections  from  the  back 
surface  of  tlic  reservoir  not  to  influence  the  motion  of  the  projectile. 

Under  these  conditions,  equation  (141  is  every  effective.  For  Ao/Ai  = 4.0 
the  maximum  error  in  the  predicted  velocity  is  about  2%  for  muzzle  velocities 
in  the  range  1.0  a®  to  2.0  ao . The  performance  for  Ao/Ai  - 1.0  is  much  the 
same,  except  near  U = 1.0  ao  where  the  error  in  U is  greater  than  2%. 
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3.  I'heory  for  iiifiiiire  reservoir  gun,  To  = 1.4 


U 

3o 


2.0  L 


1.5 


1.0 


0..5 


25.0 

4.0 


n 


® 1.0 

Data  from  .Seigel, 
figures  20, Ref. 1 


Theory  ^^y  /A,  = oo 


li 

a//- 

^ I 


bP^>/o  I 

Theory  Ao/Ai=4.0 

/© 


Q/ 
/ 


J_ 

0.5 


■Theory  Ao/Ai  = 1.0 
(equation  (14) ^ 


1.0  2.0 

po  Ai  L|  ^ 


5.0  10.0 


mao 


Figure  4.  Theory  for  infinite  reservoir  gun,  To  = 5/3 
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Similar  results  apply  for  7o  = 1.4  and  7o  = 5/3,  as  shown  in  figures  3 
;md  4.  For  7o  = 1.4,  the  muzzle  velocity  range  for  minimum  error  is  0.«  a<, 
to  1.7  while,  for  7o  = 5/3,  the  range  is  0.7  ao  to  1.5  a©.  We  conclude 
that  equation  (141  provides  a most  satisfactory  performance  when  the  driver 
gas  velocities  lie  in  the  transonic  or  the  low  to  moderate  supersonic  speed 
ranges.  Tlie  major  differences  between  figures  2,  3 and  4 may  be  removed  by 
using  the  parameters  (7oU/ao)  and  (7o  Po  Aj  Li /maS  ) . Details  are  given  in 
Appendix  III.  We  turn  now  to  the  problem  of  deciding  when  the  reservoir  is 
long  enough  for  equation  (14)  to  be  applicable. 

3.5  Effect  of  reservoir  length 

Seigel (ref. 1)  presents  some  results  for  tlie  minimum  length  of  reservoir 
necessary  to  achieve  maximum  muzzle  velocity.  Under  these  conditions,  the 
first  wave  reflection  from  the  back  surface  of  the  reservoir  reaches  the 
projectile  just  as  it  leaves  the  barrel. 

When  /V)  = Ai  , the  analytic  results  vary  only  a little  with  7oand,  as 
sho\'m  by  Sc  igel  (ref . 1)  , a satisfactory  formula  is 


_ 


ma  X i mum 


2.. 5 poA)Li/maS 


As  po  ..  0,  tlie  gas  Mach  numbers  approach  zero  and  the  minimum  reservoir 
length  approaches  infinity.  i.e.,  as  expected,  the  reservoir  volume  Vq 
for  maximum  gun  performance. 

Generalization  of  the  Ao  = A|  results  to  the  case  Ao  > Ai  presents  some 
difficulties  as  Seigel (ref . 1 ) has  given  numerical  results  for  7o  = 1.4  only. 
We  assume  that  these  results,  like  those  for  .Ao  = Ai  , arc  only  weakly 
dependent  on  7o . A simple  representation  of  reasonable  accuracy  is  then 
found  to  be 


i± 

Lo 


max i mum 


2(1-(.A, /.Ao))i  poA,  L, /m:i^  . 


(lb) 


The  value  of  lo  satisfying  equation  (IP)  is  denoted  (Lo)^.^.  Until  additional 

data  become  available,  it  is  recommended  that  equation  (16)  be  used  to  deter- 
mine when  the  reservoir  lengfli  lo  is  large  enough  for  wave  reflections  not  to 
influence  the  motion  of  tlie  projectile. 

The  problem  of  how  to  treat  cases  when  lo  is  less  th:in  (Lo)  • is  a verv 

min 

difficult  one.  In  general,  recourse  to  Seigel 's  design  charts  (1)  is  by 
far  the  most  effective  approacii.  Wlien  the  gas  velocities  arc  subsonic, 
however,  the  ijuas i -steady  theory  of  section  2 is  applicable  and  should  be 
used.  llie  limit  of  applicability  of  this  theory  is  hard  to  define  precisely. 
However,  for  gas  Mach  numbers  reaching  0.5  when  the  projectile  leaves  the 
barrel,  the  theory  is  clearly  acceptable  for  actual  guns  in  view  of  the  good 
agreement  reported  in  reference  9 between  the  theoretical  and  actual  perform- 
.ance  of  two  W.R.F..  gas  guns.  llie  maximum  muzzle  velocity  was  189  m/s,  which 
corresponds  to  a m.aximum  gas  'lacli  number  of  0.56.  Wlien  the  maximum  gas 
Mach  number  reaches  1.0  in  an  actual  gun  the  quasi-steady  theory  is  likely 
not  to  be  acceptable.  A gas  Mach  number  limit  of  0.75  is  therefore  suggested 
for  the  (pias i -steady  theory,  with  no  allowance  made  for  Lo  effects.  In  order 
to  provide  some  assistance  to  designers  at  higher  gas  Mach  numbers,  a brief 
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analysis  has  been  made  of  a representative  range  of  Seigel's  results.  As 
expected,  a simple  description  of  the  penalty  of  short  reservoirs  is  just  not 
possible.  However,  the  following  formula  may  be  of  some  use.  The  purely 
empirical  result  is 


H/lU] 


(U. 


in  in 


lo  / ( 1-0  j 


min 


1/r 


r = 4+(Ao/A,  ) 


(17) 


where  (I,))  . is  tlie  value  of  Lo  that  satisfies  equation  (lb). 

nun  IV/ 


ilie  rat  i o 


on  the  left  hand  side  of  cnuation  (17)  may  contain  errors  up  to  about  0.1. 
For  a given  Lo  , equation  (17)  sliows  that  the  velocity  loss  decreases  as  Ao 
increases . 

•An  extreme  and  unrealistic  example  of  the  limitations  of  equation  (17) 
can  be  found  by  returning  to  the  case  when  the  steady  internal  flow  model 
of  section  1.1  applies  and  the  gas  Mach  numbers  are  always  very  small. 
Equation  (9)  shows  that  Uo  constant,  when  Lq  ^ i.e.  K . In 

addition,  Uo  “ Lo  ' when  Lo  is  small,  showing  that  the  1/5  power  in  equation 
(17)  is  misleading  in  this  case.  Nevertheless,  equation  (17)  does  provide 
the  correct  limit  behaviour  in  that  LI  o as  Lo  o and  U 

^ ^ ^'-o^in- 


min 


4.  CORRECTION  FOR  -M'MOSI’HERIC  COU.NTER  PRESSURE 

So  fir,  the  analysis  lias  been  concerned  with  the  driving  pressure  force. 

Only  in  the  case  of  tiio  steady  theor>'  of  section  1.1  was  an  atmospheric  counter 
pressure  considered.  I'lie  purpose  of  the  present  section  is  to  look  at  the 
physical  basis  of  what  happens  in  front  of  the  projectile  and  then  to  derive 
suitable  correction  formulae  for  the  effects  of  counter  pressure. 

4.1  Shock  wave  and  unsteady  compression  wave  effects 

As  the  projectile  accelerates,  the  gas  in  front  of  it  is  compressed.  The 
central  problem  is  simply  how  to  determine  the  magnitude  of  this  compression 
and  how  it  varies  as  the  projectile  moves  down  the  barrel.  The  starting 
point  is  similar  to  equation  (11).  A simple  wave  can  be  set  up  to  describe 
an  unsteady  compression.  The  result  (1)  is 


P/Pi 


u 

L.  V 2 ; 

_ 

27i  (7i  -1) 


(18) 


Provided  that  the  projectile  continues  to  accelerate,  the  compression  waves 
will  coalesce  and  form  a shock  wave.  The  pressure  formula  for  a shock  wave 
i s 


P/Pi 


7i  (7i  +1) 


(19) 


1 


4 


1 5 - 
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For  Mach  numbers  u/ai  less  than  1,  the  numerical  difference  between  equations 
(18)  and  (19)  is  very  small.  Mien  u/a,  is  small,  both  equations  behave  like 

p/pi  = 1 + 7i  (u/a, ) + 0 (u’  /a,  ) . 


In  the  case  of  air,  7i  = 1.4  and  the  difference  between  equations(18)  and 
(19)  increases  up  to  .'5%  at  u/ai  = 1.  For  7t  = 1.1  and  7i  = S/7>  the  perform- 
ance is  similar.  The  conclusion  is  that  either  equation  (18)  or  equation 
(19)  could  be  used  at  subsonic  Mach  numlicrs.  Since,  however,  equation  (19) 
must  be  used  at  higher  Macli  numbers  it  is  convenient  to  use  equation  (19) 
at  all  Macti  numbers. 

Whether  a sliock  wave  forms  in  a particular  case  can,  if  desired,  be  found 
by  numerical  integration  of  the  gas  d>’namic  equations  of  motion.  In  the 
case  of  a projectile  with  constant  acceleration  (f)  , Seigel  (1)  quotes  tiie 
result  that  the  shock  wave  fonns  at  a time 


t 


shock 


2a,/  \ (7i*n  f 1 


after  the  projectile  began  to  move  and  at  a distance 


^shock 


2a]/  \ (7,  +1)  f i 


from  tlie  beginning  of  the  barrel.  At-this  time  the  proj ect i le  velocity  will 
be  f = 2ai/(7i+l)  corresponding  to  a Mach  number  of  2/(7i+l).  This 

Mach  number  is  in  the  subsonic  range,  being  at  most  1 when  7,  =1.  The 
justification  for  using  tlie  shock  wave  equation  (19)  for  calculating  the 
retarding  pressure  is  now  complete.  Only  if  formation  of  the  shock  wave 
could  be  delayed  to  liiglier  Mach  numbers  would  the  difference  between  equations 
(18)  ;ind  (19)  become  s ign  i f i c;int . 

For  tr.insonic  gas  guns,  as  distinct  from  supersonic  guns,  it  follows  that 
either  equation  (18)  or  equation  (19)  can  be  used  to  calculate  the  retarding 
pressure  acting  on  the  projectile.  Once  the  muzzle  velocity  increases  to 
supersonic  speeds,  equation  (19)  sliould  be  used. 

4.2  Steady  flow  model 

I'he  steady  flow  model  of  section  1.1  already  incorporates  a counter 
pressure  correction.  Since  the  theory  ajiplies  only  for  small  gas  Mach 
numbers,  the  additional  assuinptioti  is  made  that  tlie  muzzle  Mach  number  is 
also  small.  Hence  equation  (19)  gives  the  retarding  pressure  as  p = pi  . 

The  end  result  is  equation  (S)  or,  equivalently,  equation  (4); 


p.  A,  L,  /m  . 


This  formula  m.iy  have  wider  applicability  than  the  theory  on  which  it  is 
based. 

4.3  Method  of  Seigel  and  velocity  correction  factor 

Sei ge 1 ( re f , 1 ) has  used  a high  Mach  number  approximation  in  equation  (19) 
to  obtain 


1 + 


p/p 


7i  (7i +1  ) 


(u/a, y . 
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This  result  is  now  combined  with  a constant  driving  pressure.  i.e.  constant 

driving  acceleration.  The  simple  differential  equation  for  ^ leads  directlv 
to 


UVCUM. 


1 -il 

Pa_ 


1 - exp(-y') 


where  y = 7|  (7i  •••  llpi  A|  I,, /maf  . Ttic  form  of  equation  (20a)  is  important. 
Sul)st  i tut  ing  for  produces  a result  which  is  inferior  when  applied 

to  cases  where  the  driving  pressure  is  not  constant. 

The  right  hand  side  of  equation  (20a)  provides  a velocity  correction 
factor.  Seigel  comments  that  equation  (20a)  works  very  well  for  high 
performance  guns.  However,  for  application  to  subsonic  and  transonic  guns, 
it  should  be  noted  that  Seigel 's  result  does  not  agree  with  equation  (4)  at 
low  Mach  numbers. 

4.4  Generalized  velocity  correction  factor 

The  purpose  of  tliis  section  is  to  generalize  Seigel's  result,  equation 
(20a),  so  that  it  will  be  more  effective  for  the  lower  speed  guns  of  direct 
interest  in  the  present  report.  Combining  equation  (4)  with  equation  (20a) 
leads  to 


(HM  - ^2p,  A,  L, /m) 

Pi  =0 


l-exp(-y) 


witi'  > = 7i(7i+l)  p,  A|  l,| /ina?  . (20c) 

liquation  (20b)  is  recommended  for  design  calculations.  It  is  exact  when 
U/ai  is  small,  or  large.  IVhen  Ll/ai  is  small,  the  first  tenii  in  equation 
(20b)  dominates  and  so  the  exact  result  from  steady  internal  flow  theory, 
with  low  gas  Mach  numbers,  is  obtained.  Mien  U/ai  is  large,  the  first  term 

is  small  and  the  second  term  dominates.  Consequently,  Seigel's  high  Mach 

number  result  is  recovered  for  ll/ai  large. 

4.5  Terminal  velocity 

The  purpose  of  this  section  is  to  derive  simple  estimates  of  either  the 
terminal  or  the  maximum  muzzle  velocity  when  barrel  length  is  regarded  as 
variable.  IVhen  the  pressure  forces  depend  on  projectile  velocity  (u)  and 
not  on  distance  travelled  (x) , a terminal  velocity  is  reached  eventually. 

In  other  cases,  a maximum  velocity  is  reached  and  the  velocity  then  decreases 
owing  to  the  dependence  of  reservoir  jiressure  on  x in  the  steady  and  quasi- 
steady  internal  flow  models.  liquation  (20b)  is  not  suitable  for  calculating 
maximum  possible  muzzle  velocities  because  y depends  on  the  distance  travelled 
by  the  projectile.  We  begin  by  examining  the  steady  flow  model  of  section  1.1, 

4.5.1  Quasi-steady  internal  flow  model 

The  maximum  velocity  occurs  when  the  driving  pressure  Pj  is  equal 

to  the  retarding  pressure  p^.  Hence,  using  equations  (2)  which 

describe  the  steady  internal  flow  model. 

Pi  = Pol  1 + IA|  x/Vn)] 
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so  that  K“'  = A,  x/Vo  - (fo/pi  -1 

;uiJ  equation  (3)  gives  the  maximum  velocity'  U as 


max 


- Po^'o  j 
(7o-l)ni 


-(-ir) 


,(.7o-l)/7o 


2 PiVg  ' 

m [_  Vpi  / 


{21a) 


Ivlien  Pi  ■=  0,  the  maximum  velocity  occurs  at  x = °°  and  is 
( 2poVo  / t7o  -l)ni  1 ^ 

Equation  (21a)  covers  the  situation  where  p^  varies  with  x and  = pi . 

Another  case  where  analysis  is  possible  occurs  when  \'o  is  large  and 
ao  2^  ai  so  tliat  Pj  = jiq  , but  we  have  p^  dependent  on  velocity.  For 

this  purpose,  the  generalised  Seigel  formula  given  in  equation  (20a)  can 


be  rearranged  by  removing  (11) 
The  result  obtained  is 


pi  =0 


:md  replacing  it  by  v^po Ai  x/m  . 


= 1 2ai  (po-pi  ) / 7i  17i +npi  ] (21b) 

■Tia  A 

with  tlie  maximum  velocity  occurring  at  x = °°.  Equation  (21b)  must  be 
rejected  as  unrealistic  because  Pj  could  not  remain  constant  from 

x = 0 to  X - in  an  actual  gun.  An  alternative  approach  for  improving 
equation  (2Ia)  is  to  use  a const.mt,  mean  retarding  pressure  greater 
than  Pi  . This  is  denoteyl  by  X pi  ruul  we  assume  that  X c;in  be  determined 
by  evaluating  the  sliock  wave  equation  (19)  at  some  mean  velocity.  Hence 


X 


1 1- 

4 


U 


mean 

:ii  J 


— . + 


7i 


11 

mean 


a I 


I \ 

' 7i + 1 \ mean 

\ 4 / ai 


(223) 


,ind,  using  a root  mean  square  value  for  U from  equation  (21a)  as 
,,  mean  ' 

a first  estimate,  we  obtain 


mean 


Po  \'(l 
(7o-l)m 


X {7n-l)/7o. 

' T'  1 

Po  / 


Pi  Vo  1 
m ( 


, 1/7(1 

Po  ^ 

Pi  '' 


(22b) 


The  final  result  is 


U* 

max 


2po  Vn 

(7o  - 1 ) m 


j /'Xp,  ,(7<,-n/7o- 

Po 


(22c) 
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Equations  (22)  proviJe  an  estimate  of  the  maximum  possible  velocity 
when  the  steady  model  is  valid.  Compressibility  effects  in  the  driver 
gas  can  be  estimated  by  using  equation  (8),  derived  from  the  quasi-steady 
internal  flow  model  of  section  2.  The  velocity  IJo  calculated  from  the 
steady  model  is  replaced  by  calculated  from  equation  (22c)  to  give 


U 

\ max 


q-S 


( 


U 

max 


To 

8 


11^ 

max 


(22d) 


4.5.2  Unsteady  expansion  model 

At  transonic  iind  supersonic  driver  gas  speeds,  the  unsteady  theory 
of  section  3 is  valid.  This  is  combined  with  the  counter  pressure 
formulae  of  section  4.1  and,  since  the  counter  pressures  now  do  not 
depend  on  x but  only  on  u,  a terminal  velocity  is  reached.  The 

retarding  pressure  given  by  equation  (18)  can  be  converted  to  an 
exponential  by  assuming  that  (7i -1)  is  small,  i.e. 


^retard  " (7iu/a,  ). 

Equation  (11)  for  the  driving  pressure  when  Ao  = Ai  leads  to  a similar 
result.  Tile  case  Ao  > A|  is  given  in  e(|uation  (15),  which  leads  to 
the  exponential  approximation 


(p/pn  ) 


drive 


exp 


- 7o (u/ao ) j 


Eejuating  p . to  p,  . now  gives  the  terminal  vclocitv  II  from 

'retard  'drive  ^ max 


L. 


In  (po/pi)  = 7i  (II  ..^/a,  1 - i7o/2)  (1-(A, /Ao))  ^7o  d'  _/;><,)  . 

in  cl  A Tuii  X 


i . c.  U 


max 


(7o/2)  (1  - (A, /Aol)  + In  (po/p,  ) 
( (7()  /;io  ) * (7i  /aj  ) I 


(23) 


The  noteworthy  feature  of  this  result  is  the  logarithmic  dependence  on 
pressure  ratio.  As  previously  noted  in  discussion  of  the  unsteady 
expansion  model,  equation  (23)  Is  valid  when  the  gas  Mach  number 

is  transonic  or  supersonic.  In  addition,  the  velocity  given  by 
equation  (23)  may  not  be  realistic  for  all  reservoir  lengths.  The 

minimum  necessary  reservoir  length  has  been  examined  in  section  3.5. 

4.b  Wave  reflections  from  muzcle 

file  shock  wave  jire-.siire  formula,  equation  (19),  is  valid  up  to  the  time  at 
which  the  shock  wave  is  modified  by  pressure  reflections  from  the  muzzle. 

At  the  muzzle,  the  gas  pressure  is  dominated  by  the  atmospheric  pressure  pi  . 
The  effect  of  muzzle  reflections  can  be  allowed  for  in  a simple,  approximate 
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way  in  com]iuter  calculations  based  on  simplified  analytic  exjiressions  for 
the  driving  and  retarding  pressures.  A more  complete  treatment  requires 
full  numerical  solution  of  tlie  Inisic  gas  d\Tiamics  equations  using  the  method 
of  characteristics. 

The  first  calculation  to  be  made  is  the  position  reached  by  the  projectile 
when  the  pressure  pulse  generated  by  the  initial  movement  of  the  projectile 
down  the  barrel  has  had  time  to  travel  the  full  length  of  the  barrel  and  then 
return  to  the  projectile.  llu'  pressure  pulse  travels  at  the  speed  of  sound 


2 


Hi  . 


'Hic  second  calculation  introduces  an  approximation.  We  begin  by  noting 
that  the  pressure  near  the  gun  muzr.le  will  be  p,  , the  initial  barrel  pressure. 
It  is  now  assumed  that  a change  in  pressure  formula  must  be  introduced  in 
order  that  the  <(uasi-steady  fonnula 


P/Pi 


(7i-l)  ir 


-7i  /(7i  - Ij 


C241 


is  used  at  the  gun  muzzle.  Equation  (24)  describes  the  retarding  pressure 
when  many  muzzle  reflections  have  taken  place. 

Equation  (19)  is  used  in  full  when  the  first  muzzle  reflection  arrives  at 
the  projectile  and  equation  (24)  is  used  in  full  when  the  projectile  reaches 
the  muzzle.  In  between,  we  assume  a linear  change  in  pressure  formula 
based  on  distance  travelled  by  the  projectile  down  the  barrel. 

.S.  PEREORMAN'CE  OP  .SUBSONIC  CAS  GUNS 

We  pass  on  now  to  a review  of  the  performance  of  subsonic  guns.  The  aim  is 
to  deteimine  the  value  of  the  tlieoretical  methods  developed  in  the  previous 
sections  and,  in  particular,  to  determine  the  precision  with  which  the  muzzle 
velocities  of  actual  guns  can  be  predicted.  I'ata  arc  available  for  five  guns, 
namely  the  W.R.!;.  127  mm  gun(ref.S1  aiu!  the  W.il.l..  .S84  mm  gun(ref.9),  the  R.A.E. 

152  mm  gun(ref.l2),  tlie  N.A.E.  254  mn  gun(ref.l.S)  aiii' , as  rciiortcd  in  reference 
12,  the  CAjVTDti  152  mm  gun.  Different  methods  of  firing  projectiles  are  used. 

Ilie  two  W.R.E.  guns  employ  <iuick-.ict  ing  valves  and  are  used  primarily  for  missile 
investigations,  while  the  other  guns  emjiloy  bursting  diaphragms  .and  are  used  for 
investigating  the  effect  of  bird  impacts  on  aircraft  structures.  A wide  range 
of  sabot  designs  is  used  in  the  different  guns.  Air  or  nitrogen  is  used  as  the 
driver  gas  so  that  7o  = 7i  and  lu,  = ai  . 

No  data  are  available  for  transonic  gas  guns.  Hie  author  is  not  aware  of 
;uiy  gas  guns  operating  at  transonic  or  low  supersonic  speeds. 

5.1  I’redictcd  muzzle  velocitie-; 

Figure  5 shows  how  the  measured  projectile  speed  comjiares  with  the 
prcdicteil  speed.  In  genor.il  tenns,  the  measured  speed  is  about  10%  less 
tluui  that  predicted  by  theory.  niroughout , the  full  computer  version  of 
the  quasi -steady  theory  o'"  section  2 is  used,  including  a barrel  density 
correction,  and  wave  reflections  from  the  muzzle  are  allowed  for  as  discussed 
in  section  4.6. 
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Details  ot'  the  data  shown  in  figure  S are  given  in  table  3.  l-'igure  6 shows 
how  the  ratio  of  measured  projectile  speed  to  predicted  speed  varies  with 
reservoir  pressure.  No  clear  trend  is  ap]iarcnt.  nie  results  show  that  the 
NAi;  gun  at  its  lowest  reservoir  pressui-e  has  a measured  speed  about  20% 
greater  th:ui  the  average  tt'end.  In  view  of  tlie  fact  that  the  net  driving 
pressure  was  only  about  lui  1 f an  atiiios|ihcre  fSS  klha)  the  relatively  good 
agreement  is  [lerhaps  surjirising  and  coiifirms  tlie  value  of  the  muzzle  reflect- 
ion corrections  in  :ui  extreme  case  where  the  net  driving  pressure  is  equal  to 
about  half  the  barrel  pressure.  hlien  the  net  driving  pressure  is  about  one 
atmosphere  (105  kl’a)  , tiic  perfoi'wiance  of  the  tlieory  is  excellent. 

For  predicted  speeds  around  500  m/s,  the  results  in  figure  5 suggest  that 
the  quasi-steady  theory  is  breaking  down.  The  number  of  data  at  these  high 
subsonic  speeds  is  ratlicr  small  so  definite  conclusions  cannot  be  drawTi. 

We  can  say,  howevaM',  that  the  results  confirm  the  earlier  suggestion  that 
the  (juas i -steady  theory  can  be  used  U]i  to  a driver  gas  Mach  number  of  about 
0.75  or  so. 
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TABLE  3.  MUZZLE  VELOCITY  PREDICTIONS 


Gun 

Mass 

Pressure 

Muzzle  velocity  (m/s) 

Ratio  of  measured 
velocity  to 
predicted  velocity 

Mean  value 

(l^g) 

. 

(k  Pa  gauge) 

Theory 

Measured 

of  ratio 

2.3 

690 

111 

105 

0.95 

1.380 

163 

155 

0.95 

2070 

200 

189 

0.95 

j H-RE 

127  mm 

4.5 

690 

81 

72 

0.89 

1380 

119 

j09 

0.92 

2070 

147 

1.36 

0.93 

0.92  j 

2760 

170 

153 

0.93 

! 

1 

9.  1 

690 

56 

49 

0.  84 

i 

I 

1380 

86 

77 

0.90 

i 

1 

2070 

106 

97 

0.92 

1 

1 

1 

2760 

123 

113 

0.92 

i 

i 

1 

14.5 

345 

107 

103 

0.96 

1 

\ 

690 

163 

152 

0.93 

i 

\ 

1 

1030 

201 

184 

0,92 

i 

i 

1 

! WRE 
,384  mm 

1 

29.0 

345 

79 

75 

0.95 

i 

690 

121 

112 

0.93 

0.93 

1030 

150 

138 

0.92 

1 

1 

58.  1 

345 

57 

53 

0.93 

1 

1 

690 

87 

87 

0.94 

1 

L._ 

1030 

109 

101 

0.93 

! 

1 

i 

3.2 

55 

68 

75 

1.10 

h 1 

1 

; NAP 

103 

121 

115 

0.95 

0.98  1 

|254  mm 

175 

16  7 

155 

0.93 

; 

2 7.S 

2 i 3 

197 

0.92 

j 

1 

r 

1.8 

275 

202 

158 

0.78 

RAE 

550 

274 

229 

0,84 

0.82 

152  mm 

830 

318 

265 

0.83 

1100 

349 

287 

0,82 

CAATDC 

1.8 

760 

271 

233 

0.86 

0.83 

152  inn 

1590 

35.3 

280 

0.79 

r 
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Reservoir  pressure  (kPa  gaiiRe) 

l iRiire  ().  Influence  of  reservoir  pressure  on  performance  of  subsonic 
gas  guns 

5.2  Counter  pressure  corrections 

Section  4 was  concerned  with  the  calculation  of  atmospheric  counter 
pressure  corrections.  Tlie  purpose  of  the  present  section  is  to  examine 
the  different  corrections  discussed  and  compare  them  with  the  full  theoretical 
results  of  section  5.1.  The  difference  between  the  muzzle  velocity  without 
any  counter  pressure  and  with  counter  pressure  is  the  velocity  correction  to 
be  examined. 

An  overview  of  the  performance  of  the  generalized  velocity  correction 
factor  of  section  4.4  is  given  in  table  4.  The  maximum  error  in  the  mean 
values  of  table  4 is  25%,  and  the  minimum  is  .3%. 


TABLE  4.  MODIFIED  COUNTER  PRESSURE  CORRECTION 


1 Gun 

1 

Mean 

pressure 
(k  Pa  gauge) 

Mean 

measured 

muzzle 

velocity 

(m/s) 

Mean  value  of 
(exact  - estimated) 
velocity  correction 
(m/s) 

Mean  value  of 
(estimated/exact) 
velocity  correction 

WRE  127  mm 

1630 

115 

+ 3 

0.75 

iVRE  384  ram 

690 

111 

+ 2 

0.93 

SAL  254  mm 

150 

136 

-12 

1.12 

RAE  152  mm 

690 

235 

-10 

1.14 

CAATDC  152  mm 

1180 

257 

+ 2 

0.97 

L 
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riie  data  on  which  table  4 is  based  are  shown  in  tables  5 and  6.  Table  5 
compares  the  muzzle  velocities  for  no  counter  pressure  with  a number  of 
different  predictions  of  the  corrected  muzzle  velocities.  In  some  cases, 
the  simple  theory  for  no  counter  pressure  predicts  velocities  much  larger 
than  any  of  the  corrected  velocities.  Comparison  with  the  simple  theory 
for  the  corrected  velocity  sliows  tliat  the  performance  of  Scigel's  method 
(equation  20a')  is  encouraging  and  that  of  the  modified  method,  equation  (20b) 
iuid  table  4,  is  even  more  so.  Tlie  final  column  in  table  5 gives  the  full 
theory  with  allowance  for  the  barrel  densit)’  effect  described  in  section  2.2. 
The  difference  between  the  last  two  columns  is  small,  which  sliows  that  the 
influence  of  the  barrel  density  effect  on  the  muzzle  velocity  is  small. 

I'or  tliis  reason,  the  results  in  table  5 for  the  theory  without  counter 
pressure  were  not  re-calculated  using  the  full  theory  and  the  driver  gas 
density  correction. 


TABLE  S.  CORKICTED  MUZZLE  VELOClTIliS 


Muzzle 

■ 

ve loc i ty 

(m/s) 

Simjile  theory 

Correctec 

velocity 

('•Ull 

Mass 

I’ressure 

for  no  counter 

.Se  i ge  1 

Modi f ied 

Simji  le 

Pull 

(Kr) 

(k  I’a  gauge) 

pressure  * 

method 

method 

theory  * 

theory  ** 

2.  .3 

090 

134 

123 

112 

109 

111 

1.380 

180 

170 

162 

158 

163 

2070 

213 

204 

197 

192 

200 

4.5 

690 

96 

89 

84 

80 

81 

WRE 

1380 

130 

125 

121 

118 

119 

2070 

156 

151 

148 

144 

147 

2760 

l"7 

172 

169 

165 

170 

9.  1 

690 

(i9 

64 

60 

58 

58 

1380 

9 3 

90 

87 

85 

86 

2070 

112 

109 

107 

105 

106 

2760 

128 

125 

123 

121 

123 

14.5 

345 

150 

125 

107 

106 

107 

690 

194 

174 

161 

160 

163 

1030 

227 

208 

197 

195 

201 

WRE 

29.0 

345 

108 

91 

80 

78 

79 

690 

112 

129 

120 

1 19 

121 

10  30 

1<>7 

156 

149 

147 

150 

58.  1 

.345 

77 

66 

58 

56 

57 

690 

102 

93 

88 

86 

87 

1030 

121 

114 

109 

108 

109 

3.2 

55 

208 

108 

54 

67 

68 

NAE 

105 

233 

145 

107 

118 

121 

ZTiA  mm 

175 

26.3 

182 

151 

16.3 

167 

275 

298 

222 

194 

208 

213 

1 . 8 

275 

274 

21  1 

188 

199 

202 

RAt, 

550 

3 36 

277 

256 

267 

274 

I.SZ  mm 

830 

T’Z 

319 

299 

308 

318 

1100 

407 

350 

.329 

337 

.349 

(Wine 

1 . 8 

760 

32  5 

27*) 

258 

259 

271 

1 mill 

1590 

104 

.■^S7 

3.38 

355 

* Ixeluding  barrel  density  effect,  i.e.  gas  density  In  barrel  is  less  than 
that  in  reservoir. 

**  Including  barrel  densitv  effect. 
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TABLE  6.  COMPARISON  OF  COUNTER  PRESSURE  CORRECTIONS 


Velocity  correction  (m/s) 

Correction  ratio 

Cun 

Mass 

(ks) 

Pressure 
(k  Pa  gause) 

Seigel 

method 

Modified 

method 

Exact 

Seigel /exact 

Modi f ied/cxact 

2.3 

690 

11 

“) 

25 

0.44 

0.88 

1380 

10 

18 

*)  ^ 

0.4  5 

0.82 

2070 

9 

16 

21 

0 . 4 3 

0 . 76 

4.5 

()90 

7 

12 

16 

0.44 

0 . 75 

WRi: 

1380 

5 

0 

12 

0.42 

0.75 

127  mill 

2070 

5 

8 

12 

0.42 

0.67 

2760 

8 

12 

0.42 

0.(,7 

9.1 

690 

5 

9 

11 

0.45 

0.82 

1380 

3 

6 

8 

0.38 

0.75 

2070 

3 

5 

0.43 

0.71 

2760 

.3 

5 

7 

0 . 4 3 

0.  71 

Mean 

value 

6 

1 1 

14 

0 . 4 3 

0.75 

r 

14.5 

345 

25 

4 3 

44 

0.57 

0.98 

690 

20 

.33 

34 

0.59 

0.97 

1030 

19 

30 

32 

0.59 

0.94 

29 . 0 

345 

17 

28 

30 

0.57 

0.93 

WRE 

690 

13 

1 ■) 

2 3 

0.57 

0.96 

384  mill 

1030 

11 

18 

20 

0.55 

0.90 

58.  1 

345 

1 1 

19 

21 

0.52 

0.90 

(.90 

>) 

II 

1(. 

0.56 

0.88 

1030 

7 

12 

13 

0.54 

0.92 

Mean 

va  lue 

15 

24 

26 

0.36 

0.93 

3.2 

55 

100 

151 

141 

0.71 

1.09 

NAE 

105 

88 

126 

115 

0.77 

1.10 

254  mm 

175 

81 

112 

100 

0.81 

1 . 12 

275 

76 

104 

90 

0.84 

1.15 

Mean 

value 

86 

124 

112 

0.78 

1.12 

1.8 

275 

63 

86 

75 

0.  84 

1.15 

R/\E 

550 

59 

80 

69 

0.86 

1.16 

152  mm 

830 

58 

78 

69 

0.84 

1.13 

1100 

57 

78 

70 

0.81 

1.11 

Mean 

value 

59 

81 

71 

0.84 

1.14 

ICAATDC 

1.8 

760 

46 

67 

66 

0.70 

1.02 

[ 152  mm 

1590 

47 

66 

72 

0.65 

0.92 

' 

1 

Mean 

value 

47 

67 

69 

0.68 

0.97 

Note:  Modified  correction  ~ (3/2).  (Scigel  correction) 
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Table  6 compares  the  velocity  corrections  instead  of  tlie  muzzle  velocities 
shown  in  table  5.  The  perfomiance  of  the  Seigel  method  and  the  modified 
method  are  now  stiown  in  detail.  As  a rough  generalization,  the  modified 
method  gives  velocity  corrections  about  SO^  larger  than  docs  the  Seigel 
method.  The  very  large  velocity  correction  for  the  low  pressure  NAE  gun 
is  noteworthy. 

We  conclude  that  the  modified  velocity  correction  method  of  section  4.4, 
and  equation  (20b),  is  satisfactory  when  account  is  taken  of  the  precision 
with  which  gun  performance  can  be  estimated,  as  shown  in  figures  5 and  6. 

The  generalized  velocity  correction  factor  is  therefore  recommended  for 
design  puq>oses.  F'or  liigli  perfonnance  transonic  guns,  the  reservoir 
pressure  must  be  very  mucli  greater  than  the  barrel  pressure,  so  that  the 
velocity  correction  from  counter  pressure  is  relatively  small  and  the  magni- 
tude of  the  error  in  the  velocity  correction  must  be  a lot  less  than  the 
likely  error  in  any  theoretical  muzzle  velocity  prediction. 

5.3  Friction  and  boundary  layer  corrections  used  for  hypervelocity  guns 

Seige 1 (ref . 1)  has  compared  the  performance  of  hypervelocity  guns  with 
theoretical  predictions  aiid  so  determined  a correction.  Tl’ie  correction 
increases  as  the  driver  gas  Mach  number  increases,  and  has  been  put  down 
to  the  effects  of  projectile  friction  and  of  boundary  layer  growth  in  the 
driver  gas.  Some  careful  tests  reported  by  Seigel (ref . 1)  show  that 
variations  in  muzzle  velocity  occur  as  friction  between  the  barrel  and  the 
projectile  is  changed.  Seigel  also  reports  one  boundary  layer  calculation 
showing  that  the  effect  on  muzzle  velocity  increases  with  speed.  Both 
these  observations  are  consistent  with  the  velocity  decrements  observed  in 
the  hypersonic  gun  data.  Tlic  data  give  a velocity  decrement  of  about  2% 
wlien  TolI/;io  = 1-0  and  of  about  1%  when  7oU/ao  = 2.0.  Consequently,  when 
U/a()  is  less  than  1.0  this  effect  is  very  much  less  tlian  tlie  effects  observed 
for  subsonic  guns  in  figures  5 and  6. 

Tlie  hyjierve Jocity  gun  correction  is  tlierefore  not  helpful  in  attempting 
to  explain  the  average  velocity  loss  of  about  10%  observed  for  subsonic  gas 
guns.  However,  ballistic  gas  compressors  (6)  are  a possible  source  of 
relevant  information  and  they  are  examined  in  the  next  section. 

5.4  Corrections  used  for  ballistic  gas  compressors 

In  ballistic  gas  compressors,  a moving  piston  is  used  to  compress  g;is  in 
a sealed  tube.  The  kinetic  energy  of  the  piston  is  utilised  to  heat  the 
gas.  Alkidas,  I’lett  and  Suinmerfield(ref .6)  liave  identified  a number  of 
different  effects  that  must  be  included  in  any  method  of  predicting  the 
performance  of  ballistic  compressors.  The  import.int  effects  are 

(a)  real  gas  effects  (e()uation  of  statel  , 

(b)  valve  losses  at  entrance  to  barrel, 

(c)  gas  leakage  between  piston  and  barrel, 

(d)  generation  of  shock  waves  in  the  compressed  gas, 

(e)  heat  losses,  and 

(f)  friction  between  piston  and  barrel. 

Items  (c)  ;ind  (f)  are  the  only  items  relevant  to  transonic  gas  guns.  .All 
the  other  items  arise  either  because  of  the  very  high  pressures  and  temp- 
eratures generated  in  the  ballistic  compressor  or  because  of  special  features 
in  the  ballistic  compressor. 

The  major  causes  of  the  discrepancy  observed  between  theory  and  experi- 
ment in  figures  5 and  6 are  taken  to  be 

(i)  gas  leakage  between  barrel  and  sabot,  and 

(ii)  friction  between  barrel  and  sabot. 
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lil't’ect  (i)  is  always  present  to  some  extent  and  a method  for  calculating  it 
is  yivon  in  reference  6.  The  method  is  based  on  incompressible,  viscous  flow 
and  requires  as  input  data  the  size  of  the  gap  between  the  sabot  and  barrel, 
and  the  pressure  difference  across  the  length  of  the  sabot.  Effect  (ii) 
cjuuiot  be  calculated  and,  if  necessary,  may  be  best  determined  empiricall>'  as 
it  will  probably  vary  a lot  from  sabot  to  sabot  and  from  gun  to  gun.  Some 

tests  on  the  W.R.E.  .384  mm  gun  show  tltat  effect  (ii)  was  not  important  for 
that  gun.  However,  this  gun  did  not  have  a machined  barrel  so  that  effect 
(i)  may  have  been  significant. 

5..S  Empirical  assessment  of  subsonic  gun  performance 

Figures  5 and  6 show  that  the  actual  projectile  speed  is,  in  general  terms, 
about  10°ii  less  than  predicted.  There  seems  to  be  no  simple  way  of  improving 
the  theoretical  prediction  methods.  The  discussion  in  section  5.3  ;ind  5.4 
has  led  only  to  the  i|ualitative  conclusion  that  gas  leak.age  .and  friction 
between  tile  s.abot  and  b.irrel  are  the  main  causes  of  the  10%  or  so  loss  in 
performance.  Conseiiuent  ly , for  design  purposes,  we  m.ake  the  empirical 
assumption  that  all  theoretical  velocities  should  be  reduced  by  10%  as  a 
final  step  in  predicting  projectile  speeds. 


6.  SUMMARY  01  DF.SION  PROCEDURE 

When  the  Mach  numbei'  in  the  driver  gas  is  subsonic,  U'ss  than  0.75  or  so.  the 
quas  i - s t I'.idv  theory  .applies.  Die  c.ilculation  steps  are  the  muzzli'  velocity  II 
when  theie  is  lu)  gas  in  the  Inirrel  i n i t i . 1 1 I y , a counter  pressure  correction  for 
gas  in  the  b.irrel  and,  finally,  the  empirical  10%  correction  of  section  5.  I'he 
simple  analytical  procedure  is  summarized  in  table  7 below. 


TABLE  7.  SUBSONIC  UUS  DESION,  U/;io  < 1 


-Step 

number 

Vo  loc i t y 

Formula 

(aimment 

1 

Uo 

/ "’ll  \ I k'  k'4-  1 ” ^ ^ 

( lJo/3o)  = -{’'•''‘-i.  . 1-/  * 

^ ^ ma<^,  (7n-l)  L \ K J 

eipiat  i ons 
(3), (9) 

1 J 

1 

1 

(U) 

Pi  =o 

y 1 1/  0^  ~ ('  ^)  / ■>()  1*1  1 ~ ^ 1 Ifi  /*lo  ) 1 

\ Ti  1 1 

eiiuation  (8) 

3 

U 

1''  = =o  * f^p.  A,  L, /mjj. 

equat i on 
(20b) 

where  y = Ti  (7i  <-npi  A|  L|  /maj 

equat ion 
(20c) 

4 

1- 

11 

est imate 

U _ ^ = 0.9  11 

estimate 

empi r i ca 1 
correction , 
section  5.5 

L 
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When  the  Mach  number  in  the  driver  gas  reaches  transonic  or  supersonic  speeds 
a different  procedure  is  used  to  calculate  the  muzzle  velocity  U when  no  gas  is 
in  the  barrel.  The  entire  procedure  is  summarized  in  table  8 in  a form  suitable 
for  immediate  use  by  a designer. 


rABLb  8.  TRANSONIC  OR  SUPERSONIC  CUN  DESIGN,  U/ao  ^1 


EX/\MPLE; 


In  order  to  illustrate  the  use  of  the  analytic  design  procedures  summarized 
in  section  6,  the  performance  of  a liypothet ical  transonic  gas  gun  has  been 
calculated.  Figure  7 and  table  9 compare  the  performance  of  h)'drogcn , helium 

and  nitrogen  as  driver  gases.  Nitrogen  is  markedly  inferior  to  hydrogen  or 
helium  because  the  driver  gas  Mach  numbers  arc  transonic  and,  as  a direct 
consequence,  there  is  a substantial  drop  in  driving  pressure  as  the  projectile 
accelerates.  On  figure  7,  assuming  a design  Mach  number  of  1 . X2 , a nitrogen 
gun  witli  a 27t  m barrel  has  the  same  pcrform.ince  as  a helium  gun  with  a barrel 
10  m in  lengtli.  For  the  same  performance  a hydrogen  gun  would  have  a slightly 
shorter  barrel,  8.8  m long.  For  the  2X  m barrel  with  nitrogen,  the  reservoir 
length  (and  volume)  have  been  increased  slightly  in  order  to  satisfy  the  length 
requirement  of  table  8. 
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Table  9 compares  the  mass  (H)  of  gas  in  the  reservoir  with  the  projectile 
mass  (m) . For  the  cases  considered  in  the  table,  the  ratio  G/m  is  always  a 
little  less  than  unity  for  helium  and  hydrogen,  but  is  markedly  greater  than 
unity  for  nitrogen.  Hence,  if  using  nitrogen  as  a driver  gas,  transonic  muzzle 
velocities  demand  that  the  mass  of  driver  gas  be  greater  than  the  mass  of  the 
project i le . 


Muzzle 

Mach 

number 


1.78, 


1 - .8 


1.28 


1.0 


0.  :’8 


/ /" 

Hydrogen  p / 

// 

// 

' ' Hel i urn 


/ 


7 


/ 


X 


/ 

/ 

'^Ni trogen 


(Vo=0.89in'  j^ 


Po  - 1.8  MPa 
A,  =0.02  in’ 
m = 1 8 k g 
V„  =0.8  m' 
A,i  =0.1  m’ 


0 10  20 
Barrel  length  fm) 


.80 


Figure  7. 


Fxample  of  transonic  gun  design 
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TABLE  9.  CUN  DESIGN  EXAMPLE 


Basic  gun 
parameters 

Driver 

gas 

Driving 

par.'uneter 

/ po  Ai  Li  'v 

Retarding 

parameter 

Ti  (Ti  l)pi  A)  Li 

Muzzle 

velocity* 

(m/s) 

Muzzle 

mach 

number* 

k mao 

ma| 

i 

1 

1 Reservoir  press  ire, 

I |\)  = 15  MPa 

Barrel  area. 

Nit  rogen 
an  = 350  m/s 
7o=1.4 
(;/m=6.72 

1.033 

0.0  ?i8 

1.08 

1 

j 

! A,  =0.02  m^ 

1 

1 

Barrel  length, 

L,  = 10  m 

Helium 
ao=1000  m/s 
To =5/3 

G/m=0. 83 

0.200 

0.038 

456 

1.32  1 

i 

1 lotal  mass, 
m=15  kg 

j 

: Reservoir  volume, 
V„=0.5  ra-' 

Reservoir  area, 

/\o=0.  1 in’ 

Hydrogen 
ao  = 1300  m/s 
To  = 1 • 4 

G/m=0. 41 

0.  118 

0.038 

482 

L 

1.40  ^ 

; 

Hydrogen 

L)  =5  m 
G/m=0.41 

0.0  S9 

0.019 

367 

1.06  ! 

1 

; 

Air  pressure, 
p,  =0.  1 MPa 

Air  speed  of 
sound 

ai =345  m/s 

He  1 i urn 

Li  = 5 m 
G/m=0 . 83 

0.  100 

0.019 

1 

1.0,^  i 

1 

1 

Helium 

L|=15  m 
G/m=0. 83 

o 

0.057 

513 

1.49 

Air  ratio  of 
specific  heats, 

Ti  = 1.4 

1 _ 

Nitrogen 

Li  = 15  m 
C/m=5. 72 

2 . 450 

0.057 

413 

1.20 

Nitrogen 
Li=23  m** 
r./m=5.72 

3.756 

1 

0.087 

456 

1 

' 

rsj 

t-o  , 

r-H  1 

! 

* Emjiirical  loss  factor  of  10°i)  is  included. 

“ l or  nitrogen,  Li  =2.'S  m is  equivalent  to  either  po  = .^2  MPa  or  m=7.0  kg.  The 
case  Li =23  m requires  Vo=0.59  m’  so  that  the  reservoir  will  bo  long  enough 
(5.9  m)  to  stop  reflections  from  reaching  the  projectile. 


figure  8 and  table  10  show  how  gun  perfoiinajice  varies  with  reservoir  pressure, 
nie  extreme  cases  of  hydrogen  ;uid  nitrogen  are  compared  and  show  that  there  is 
a marked  increase  in  the  performance  advantage  of  hydrogen  as  the  reservoir 
pressure  increases.  Tlie  performance  of  helium  would  fall  a little  below  that 
of  hydrogen,  being  very  much  better  than  that  of  nitrogen.  Figure  8 shows  that 
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the  high  Mach  numbers  in  tlie  case  of  nitrogen  do  indeed  have  a dramatic  effect 
on  muzzle  velocity.  As  ttie  pressure  drops,  the  performance  of  nitrogen  and 
h\ drogen  becomes  much  less  different.  In  the  limit  of  low  pressures,  the 
nitrogen  and  liydrogen  guns  must  have  identical  muzzle  velocities  as  compressibility 
effects  in  tlie  driver  gas  are  no  longer  significant. 


Reservoir  pressure (MPa) 


figure  8.  Example  showing  performance  variation  with  reservoir  pressure 
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TABLE  10.  EFFECT  OF  RESERVOIR  F’RESSURE 


Driver 

gas 

Reservoir 

pressure 

(Ml’aj 

I'r  iving 
par;unetcr 

PoA|  Lr^ 

Retarding 

parameter 

7i  (7i  1 j pi  Ai  Li 

Muzzle 
velocity  * 
(m/s) 

Muzzle 

mach 

number  * 
(H/a,  ) 

; j 1 

. mao  '' 

mat 

10 

0.079 

.397 

1.15 

15 

0.  1 18 

482 

1.40 

Hydrogen 

20 

0.  158 

0.038 

551 

1.60 

30 

0.237 

660 

1.91 

40 

0.  31b 

744 

2.16 

10 

1.088 

326 

0.94 

15 

1.633 

371 

1.08 

Nitrogen 

20 

2.  177 

0.038 

403 

1.17 

30 

3.265 

449 

1..30 

32 

3.4  83 

456 

1..32 

40 

4.354 

481 

1.40 

‘Empirical  loss  factor  of  10°4  is  included 

Note:  Cun  parameters  are  taken  from  gun  design  example,  with  Li  = 10  m 
and  Vo  = 0.5  m throughout. 


S.  CONCLUSIONS 

The  primary  aim  of  this  report  has  been  to  derive  simple  analytic  design 
formulae  in  forms  suitable  for  the  use  of  gas  gun  designers.  The  main  conclus- 
ions are  listed  below. 

(1)  The  case  of  subsonic  flow  by  the  driver  gas  can  be  treated  analytically. 

A simple  analytic  formula  based  on  a quasi-steady  theory  has  been  derived. 
Muzzle  velocity  results  from  gas  guns  show  that  the  theory  is  applicable 
up  to  a gas  Mach  number  of  0.75  or  so. 

(2)  Mien  the  driver  gas  moves  at  transonic  or  supersonic  speeds  a simple 
design  formula  can  be  derived.  Hie  range  of  usefulness  extends  up  to  a 
Mach  number  of  2.  Tlie  basis  of  the  formula  is  an  unsteady  simple  wave 
expansion.  Tlie  effect  of  an  area  change  at  the  barrel/reservoir  junction 
is  included  in  the  general  formula. 

(3)  There  is  a minimum  reservoir  length  for  the  validity  of  the  unsteady 
expansion  wave  model.  A simple,  empirical  expression  for  the  minimum 
length  has  been  derived.  Tlie  unsteady  theory  applies  when  the  reservoir 
is  sufficiently  long  to  prevent  the  expansion  wave  reflected  by  the  back 
of  the  reservoir  from  reaching  the  projectile  before  it  leaves  the  gun 
barrel . 

(4)  It  is  very  difficult  to  estimate  the  muzzle  velocity  when  the  reservoir  is 
too  short  for  the  simple  unsteady  tlieory  to  be  valid.  A simple,  analytic 
formula  of  wide  applicability  cannot  be  found.  For  this  reason,  a formula 
of  limited  value  has  been  derived  empirically  to  give  some  indication  of 
the  loss  in  performance  with  short  reservoirs. 
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(51  nie  effects  of  counter  pressure  can  be  treated  analytically.  A simple 
correction  formula  has  been  deducctl.  The  coirection  is  recommended  for 
.yoTieral  use  ;uid  has  been  derived  so  that  it  a^;r' e noth  with  Seij;cl's 
correction  formula  wlien  the  mutzle  velocitv  is  1.,  .e  and  with  the  simple 
steady  theory  when  the  Mach  number  is  low. 

(6)  Estimates  can  be  made  of  the  j;reatest  velocity  possible  in  a jiarticular 
i;iui , when  barrel  lenjtth  is  treated  .is  :i  variable  ipi.int  i t\'.  A maximum 
velocitv  i'ormul.i  lias  been  dcriveil  b\'  allowing  li.irrel  lenyth  to  v.ary  in 
the  ipias  i -steady  theory.  Simil.irli',  .a  terminal  velociti’  formul.i  has 
I'cen  derived  for  higher  speeds  where  lioth  the  driviii).;  and  ret.ardiriK 
pressure  are  dependen  'n  velocity  only. 

(7)  The  performance  of  actual  gas  guns  is  less  thati  that  predicted  b\'  theory. 

An  analysis  has  been  made  of  tlie  performance  of  subsonic  gas  guns  in  order 
to  obtain  guidelines  for  the  likel>'  errors  in  theoretical  muzzle  velocity 
predictions.  Overall,  the  velocity  is  about  lO^o  loss  than  predicted. 

The  loss  in  performance  is  attributed  to  gas  leakage  between  the  sabot 
and  the  barrel,  ;uul  to  friction  between  the  barrel  and  the  sabot. 

(8)  hither  hydrogen  or  helium  is  superior  to  nitrogen  as  a driver  gas  for 

transonic  guns.  A design  example  sliows  that,  for  similar  muzzle  velocities, 
tlie  reservoir  pressure  in  the  case  of  nitrogen  is  about  twice  that  required 
for  liydrogcn  or  helium.  Because  the  speed  of  sound  in  hydrogen  is  greater 

tiuui  that  in  heliimi,  hydrogen  .ilw.ivs  gives  slightly  better  performance  than 
he  1 i urn . 
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N'O  TA  r 1 ON 

A 

spocd  of  souiul 

(a) 

son  1 c 

value  of  'a'  wliere  the  yas  velocity  is  equal  to  'a' 

A 

area  of  cross-section 

C 

const  ;uit 

u 

Ji  luneter 

f 

acceleration  of  projectile,  when  constant 

c 

K 

mass  of  i^as  in  reservoir 

1 

= V„  /A,  L, 

i !. 

■ 

leiiKtli  of  reservoir  or  barrel 

(Ih)1  . 

mill 

minimum  leiu',th  of  reservoir  to  avoid  wave  reflections  at  projectile  | 

m 

mass  of  projectile,  includiiij^  sabot 

I’ 

jiressure  of  gas 

I’J 

driving  jiressure 

retarding  pressure 

r 

- 4 »■  (Ai)/Ai  ) 

t 

time,  mcasureil  from  when  projectile  first  begins  to  move 

1 

r ^'iliock 

time  t at  which  shock  wave  is  first  formed 

1' 

temperature  of  g.as 

11 

velocity  of  proiectile  in  barrel 

II 

miiccle  ve  loci  tv 

lln 

value  of  U when  ('H/a,,  ) 0 (steadv  theory) 

u 

ost  imiite 

estimated  value  of  H for  a gas  run 

11 

max 

maximum  value  of  U when  L|  varies 

11 

mean 

mean  value  of  U,  equation  (22b) 

""p.-O 

value  of  11  when  pi -0 

value  of  U when  U = (Lo)  • 

nun 

v„ 

volume  of  reservoir 

X 

distance  travelled  by  projectile 

■''■-lioek 

position  ot'  shock  wave  when  first  formed 

\ 

retarding  parameter  7i  (7i +l)pi  A)  I,| /mat 

= ( /(ToHl/f  - 1)  ll  - (A,  AVj) 

L 
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7 ratio  of  specific  heats 

^ ratio  of  mean  retarding  pressure  to  initial  retarding  pressure  pi 

P density  of  gas 

Subscripts 

0 initial  reservoir  conditions 

1 initial  barrel  conditions 

s gas  stagnation  conditions 

q-s  ipias i -steady  model 
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Ai’I’liNlUX 


THE  PIDDUCK-KENT  Sl’ECrAL  SOLUTION 


When  the  mass  (0)  of  gas  in  the  reservoir  is  very  small  compared  with  the 
projectile  mass  (m) , the  "I’iddiick-Kent  Special  Solution",  often  used  in  internal 
ballistics  and  described  in  reference  1 for  example,  predicts  the  same  muzzle 
velocity  as  the  steady  theory  of  equation  (3).  Tliis  is  not  surprising  because 
the  "I’idduck-Kent  Special  Solution"  is  a most  useful  approximation  of  wide 
applicability  (1).  However,  in  order  to  obtain  an  analytic  solution,  it  does 
assume  that  there  is  a pressure  gradient  of  a particular  form  and  this  assumption 
could  lead  to  a different  quasi-steady  theory.  The  muzzle  velocity  correction 
for  higher  speeds  can  be  found  from  the  series  expansion  (1)  in  powers  of  (G/m) . 
Tlie  result  for  small  values  of  (G/mj  is  a velocity  correction  factor  equal  to 
1 - I (37o-l)  /12'yo  i (G/m) , which  is  different  from  the  quasi-steady  correction 
of  equation  (8).  In  both  cases,  the  small  correcting  term  is  proportional  to 
Uo/ao.  From  equation  (9), 


Uo  / ao 


7o  (7o-l) 


-(7o-l) 


and  the  velocity  correction  factor  of  equation  (8)  can  be  written  as 


(G/m)  ) K+ 1 '\  - ('I'o  - ^ ^ 

4(70-1)  ' 'v  K ; 


Therefore  the  quasi-steady  correction  depends  on  K as  well  as  7o  and  (G/m)  whereas 
the  correction  from  the  "Pidduck-Kent  Special  Solution"  depends  solely  on  7o 
(G/m).  The  coefficients  of  (G/m)  arc  respectively 


4(7o-l) 


1 

V K 


- (7o-l) 


(37o-1) 
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APPliNDIX  II 

SKU'.EL  PRHSSUUP.  POUMllLA  POK  An  > A, 

llie  barrel  entry  sonic  approximation  is  descrit>ei.i  in  section  3.3.  I.quation 
(15)  gives 


P/Po 


1 + Z 


7o_-i)  u_ 

^ Z / 


Z7u  / (7o  - 1 ) 


(II.  n 


where 


j^\/f7o+l)/2 


I - (A, /An) 


Seigel  (ref . 1)  lias  introduced  an  iinjiroved  fonnula  h\-  tal'.ing  equation  (II.l)  to 
apply  only  for  70^/110  > 1.5.  Mu'n  7()U/;jn  < 1.5,  Z is  replaced  by  (7o  ti/1 . 5;:o  ) . Z 
which  leads  to  the  correct  behaviour  in  the  limit  as  u 0.  Seigel's  improved 
formula  leads  to  a muzzle  velocit>'  formula  more  complex  th;\n  equation  (14). 
However,  in  cases  where  tlie  distance  travelled  for  a given  muzzle  velocity  is 
required,  Seigel's  improved  formula  is  st  ra  i ght  foi'wanl  to  use  in  the  equation  of 
motion  (6).  'Die  inversion  to  g.ive  velocity  explicitly  is  not  so  straightforward 
and  must  be  restricted  to  transonic  and  su|iersonic  speeds  if  a logarithmic 
formula  is  desired.  In  view  of  the  good  jicrformance  of  ei^uation  (14)  at  tr;ui- 
sonic  and  supersonic  speeds,  as  described  in  section  3.4,  Seigel's  more  complex 
pressure  formula  has  not  been  used.  nie  simplicity  of  equation  (14)  is  worth- 
while preserving. 

The  case  7oU/a<)  1.5  will  now  he  exiuiiined  in  more  detail.  T)ie  factor  Z 

in  equation  (II.  1)  is  multiplied  by  (7o  u/ 1 . 5a()  ) . I'he  resulting  pressure  formula 
therefore  looks  like  equation  (11)  with  ])(,  unchanged  and  ao  replaced  by 


>'U) 


1 - 


^oZ 

1.5(7o-1)._ 


llie  a,)  multiplying  factor  becomes 


A 

1 1 - (1-  I bv  introducing  an  expansion  based  on  (7o-l)  being  small, 

I 3 / V)  ] 

Hence  equation  (13)  for  the  muzzle  velocity  becomes 


U/ao 


_ 1 

27, , I 1 - (7o73)  (T-  (A,  /,vr) ) I 


. In 


7„  , Ap 


Po  Ai  Li  /nao 


1.0  <7oH/ao  < 1.5. 


IVlien  7oH/:i<>  ^ 1.5,  even  more  comtilexity  will  be  introduced  in  the  formula.  Thus 
Seigel's  improved  pressure  formula  is  best  used  to  calculate  disttince  travelled 
for  a prescribed  muzzle  velocity.  nu;  result  for  7oll/a<,  >1.5  can  be  shown  to  be 


ll/a,. 


_L  ■'  Al\ 

47o  ^ Ao  / 


+ . In  [ i;+((i7?,  p„A,  L, /ma/,  )l  , 7o  H/tto^  1 . 5 , 
-7,1 


where  li 


exp  I 3-(l/2)  (1-(A, /A„))|  -J  l-f7o/3)(l-(A, /A«))  pexpf3-7n(Ai/Ao))l. 
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APPENDIX  in 

SIMILARITY  PARAMETERS  FOR  UNSTEADY  EXPANSION  MODEL 

Tile  pressure  formula  in  an  unsteady  expansion  is  given  by  equation  (11), 
namely 


p/po  = 

~ 1 ^ 

\ 2 / ao_ 

27o/(7o-l) 


Treating  (7o-l)  as  small,  this  formula  becomes 


p/po  = exp  (-7ou/ao)  . 


CIII.l) 


Seigel (ref . 1)  comments  that  equation  (111.1)  is  a useful  simplification.  Further- 
more, it  shows  that  7ou/;io  is  likely  to  he  an  important  similarity  parameter. 

We  can  now  combine  equation  (III.l)  with  tlie  equation  of  motion  (6)  to  obtain 

m u -T-  = A,p„  exp  (-7ou/afl ) . 


(7ou/ao)  exp  (-7ou/ao) 

d(x/l|)  mae 


Hence  we  expect  tliat  the  only  similarity  parameters  of  importance  in  determining 
mu7,zle  velocity  are  likely  to  be  7nU/a<)  and  7oPoAi  Lj /mao  . Tito  area  ratio  Ao/A|  , 
or  the  diameter  ratio  Do/D(  , completes  the  set  of  parameters. 

Hie  usefulness  of  the  similarity  parameters  7oU/ao  and  7opoAiL)/mao  can  best 
be  determined  by  re-plotting  the  data  given  in  figures  2,  3 and  4.  Tlie  results 
are  given  in  figure  11  I.  1 and  show  that  the  similarity  parameters  are  remarkably 
effect i ve. 

Tlie  predictions  of  the  simple  logarithmic  formula,  equation  (14),  arc  shown 
in  figure  III.l  together  with  data  fur  Ao/A|  ;=  1,25  taken  from  Seigel's  exact 
resul  t s (ref . 1 ) . flic  slope  of  the  str.iight  lines  representing  equation  (14)  is 
a good  average  representation  of  the  data.  For  Ao/A|  - 1,  the  slope  is  a little 
larger  than  optimum  while,  for  Ao/A|  - 25,  the  slope  is  a little  less  than  the 
optimum.  IMgure  III.l  therefore  confiniis  the  usefulness  of  equation  (14)  as  a 
sim|)le,  eftective  design  formula. 
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liS'.iiro  II  1.1  Use  ot  similcirity  parameters  1 nr  inlinite  reservoir  (jiin , 

■y,.  = l.I,  l.‘I,  .S/A 

Another  feature  of  figure  III.!  is  tlie  very  gooil  agreement  hetween  tite 
simple  logarithmic  formula  and  Seigel's  dat.i  for  A„/A,  = 25,  with  y.U/ao  < 2. 
Otis  is  significant,  because  it  covers  the  region  of  most  interest  for  guns 
operating  at  transonic  speeds.  Mten  A<,/Ai  = 1,  the  logarithmic  formula  is 
very  effective,  but  not  quite  so  good  as  in  the  Ao/Ai  = 25  case. 
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